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Abstract 

Observational limits on the high-energy neutrino background have been used to place general 
constraints on dark matter that annihilates only into standard model particles. Dark matter 
particles that annihilate into neutrinos will also inevitably branch into electromagnetic final states 
through higher-order tree and loop diagrams that give rise to charged leptons, and these charged 
particles can transfer their energy into photons via synchrotron radiation or inverse Compton 
scattering. In the context of effective field theory, we calculate the loop-induced branching ratio 
to charged leptons and show that it is generally quite large, typically > 1%, when the scale 



of the dark matter mass exceeds the electroweak scale, My/. For a branching fraction 



> 



the synchrotron radiation bounds on dark matter annihilation are currently stronger than the 
corresponding neutrino bounds in the interesting mass range from 100 GeV to 1 TeV. For dark 
matter masses below My/, our work provides a plausible framework for the construction of a model 
for "neutrinos only" dark matter annihilations. 
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While dark matter (DM) accounts for 25% of the total energy density in the universe (a 
classic review is [1]), direct detection remains elusive (see, e.g., Ref. [2] for a recent discus- 
sion). If the dark matter is a thermal relic, one would expect a non-neglible annihilation cross 
section, and one might hope to either detect the products of such annihilations occurring 
today, or to bound the properties of such particles by the non-observation of annihilation 
products. 

Recently, Beacom, Bell and Mack [3] made an interesting argument regarding such limits. 
They noted that if a dark matter particle (x) couples only to Standard Model (SM) particles, 
then the most general limits can be derived by assuming annihilation into neutrinos, as 
these are the most difficult SM annihilation products to detect. Rather surprisingly, current 
observations place nontrivial bounds on such annihilations (see also Refs. [4, 5]). Mack et 
al. [6] compared these neutrino annihilation limits with corresponding limits on gamma-ray- 
producing annihilations. Their results indicate that, for an illustrative branching ratio to 
photons of 1CT 4 , the neutrino bounds generally provide tighter constraints for large dark 
matter masses (M x > 100 MeV), while the photon bounds dominate at smaller masses. 
If the branching ratio into photons is taken to be larger, then the mass range over which 
the neutrino bounds dominate shrinks accordingly. Further, with the GLAST satellite [7] 
expected to launch in 2008, the photon limits will shortly be considerably tightened (or dark 
matter annihilation will be detected!). 

Thus, it is very interesting to determine what constitutes a "reasonable" minimum 
branching ratio into electrically-charged particles, for these in turn generate photons. We 
will assume a "neutrinos only" final state at tree level, and then calculate the branching into 
charged leptons that results from a higher-order box diagram. Note that a calculation similar 
in spirit was undertaken by Kachelriess and Serpico [8], who examined the electromagnetic 
mode from electroweak bremsstrahlung of real W and Z particles. In quantitative detail [9] , 

M 2 a M 2 

it is found that the branching fraction of this process is = 2.1 x 10 x for 

M x > M w such that the electroweak bosons are produced on-shell; we have taken a = 1/128 
for the numerical value, as is appropriate at the weak scale. The processes we discuss here 
are physically distinct from brehmsstrahlung and offer a potentially more significant con- 
tribution to the electromagnetic branching ratio. In this work, we obtain a rate ratio from 
loop graphs growing as ( 8 "^2 j times a logarithmic factor. 

It is non-trivial to embed a tree-level "neutrinos only" final state model into a field theory. 
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The special status afforded the neutrino breaks the SU(2) invariance of the weak interaction. 
As a consequence, the embedded theory we construct must be viewed as incomplete, as an 
effective theory at best, valid up to some energy scale A. The effective theory will also 
be non-renormalizable, and therefore unstable against radiative corrections. This presents 
another indication that the neutrinos-only model exists only for a finely-tuned effective 
theory below a scale A, or for a UV-completed theory with a rich spectrum at A ~ Mb, 
where Mb is the mass scale of the interaction connecting the DM sector to neutrinos. The 
rich-spectrum possibility seems contrived and unlikely, but if true, it predicts much discovery 
not far beyond the energy reach of today's accelerators. We will assume that the neutrinos- 
only model is simple in a region of validity up to A 2 ^> M B . We will present the consequences 
of such a model. 

Assuming the dark matter to be fermions, we must at a minimum introduce a new scalar 
or vector boson particle/field B, with mass Mb, to mediate the dark matter annihilation to 
a neutrino-antineutrino pair. For the process x + X ~~ *■ v u (here we are using the generic v to 
include all three neutrino flavors), one can imagine two types of vertices and diagrams: either 
an s-channel tree graph (I) with vertices mediated by C xx b and L v vb, or a t-channel tree 
graph (II) with both vertices mediated by C x b v - A new quantum number is needed to elevate 
the xxpv operator to a privileged status, but we do not indulge in detailed model building 
here. Rather, we focus on the Standard Model (SM) loop-corrections to the diagrams (I) 
and (II). 




(I) s-channel B-exchange (II) t-channel B-exchange 



There are order (g 2 ) radiative corrections to the s-channel process which produce two- 
particle charged leptonic final states Z + Z . Since we are comparing to the two-particle 
final states Xw^W> the number of flavor generations will cancel from our results. For 
example, the ratio is unchanged whether the new interaction is assumed to produce all 
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three neutrino flavors, with concomitant induced production e + e~, t + t~, or whether 

the new interaction is assumed to produce, say, just v T v T with concomitant production of 
just t + t~. Corrections to the two-particle s-channel process are (a) s-channel Z-exchange, 
and (b) t-channel W-exchange. (These SM processes necessarily violate conservation of any 
charge carried by the B particle.) 




(a) oblique correction (b) induced Bl + l vertex 



The "oblique" correction in (a) induces B-Z mixing. One may choose a counterterm to 
eliminate the mixing at one scale fi. Although not natural, one can in principle choose the 
scale fi to eliminate the mixing at the threshold fi = 2M X ; then, for non-relativistic DM, 
there is almost no production of from non-relativistic \ annihilation. Similarly, the 
Bl + l~ vertex correction in (b) can be canceled by a counterterm at one scale, e.g. fi = 2M X . 

Radiative corrections to the t-channel process which produce an final state are an 
s-channel Z-exchange and a t-channel ^-exchange. 




(c) induced Z\X vertex (d) induced xxvv box graph 



In graph (c), the radiatively induced XX.Z vertex can be canceled by a counterterm at 
one value of fi, say fi = 2M X . In contrast, the four-fermion vertex induced in graph (d) 
is of dimension six and has no counterterm. Thus, it must be rendered finite by direct 
calculation. The superficial degree of divergence in graph (a) is quadratic, (b) is quadratic 
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(in unitary gauge), (c) is logarithmic for scalar B and quadratic for vector B (in unitary 
gauge) and (d) is logarithmic for scalar B and quadratic for vector B (in unitary gauge). 
We now focus our attention on graph (d). We do so for two reasons. The first reason is that 
there is no counterterm for this graph to suppress the operator, and the second reason is 
that graph (d) is the least divergent of the graphs (along with (c)) and therefore represents 
a conservative window to our study of the stability of a "neutrinos only" model in the face 
of radiative corrections. 

For our calculation, we make the simplest choice for the spin-parity of the £?-meson, 
namely that it is a scalar particle. In fact, we show in Appendix I that the divergent 
contribution to the amplitude of the box graph factorizes into a divergent factor independent 
of the 5-meson's spin-parity, times a factor proportional to the tree-level amplitude. This 
latter factor depends on the £>-meson's spin-parity. However, this latter factor cancels in 
the ratio of box to tree amplitudes, which is what we investigate. Thus, our results are true 
for any spin-parity of the 5-meson. 

We render the box graph (d) finite by use of a Pauli-Villars (PV) regulator to compensate 
the B particle propagator. The Pauli-Villars cutoff A then becomes the maximum scale for a 
credible effective theory. We work in the unitary gauge, since this gauge, and only this gauge, 
includes only physical particles throughout the calculation. Thus, the only new parameter 
introduced is the PV regulator A, which now takes on the significance as the mass-scale 
above which additional new particles enter. Some details of our calculation in the unitary 
gauge are presented in Appendix (I). In Appendix (II) we calculate the same box graph in 
the general i?^-gauge, and compare the results to that from the unitary gauge calculation. 

The amplitude for graph (d) depends on four dimensionful parameters, M x , M B , Mw, 
and A. Therefore, the ratio of rates (independent of the number of flavor generations) 

n _ {v<XX^l + l-)) (1) 
(va(xx -> v v )) 

depends on three independent ratios, which we may take to be M x /M w , M B /M X , and A/M B . 
My/ is known, of course. For the case of a thermal relic (only), M x and M B are constrained 
by the requirement that the x dark matter decouple from early universe thermal equilibrium 
such that VLouh 2 = 0.1, where h is the Hubble parameter in units of 100 km sec -1 Mpc -1 . 
We will investigate this constraint below, in Eqs. (2-5). 

The PV mass A is a priori a free parameter. Inherent in the PV regularization scheme 
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is vanishing of the electromagnetic amplitude at A = M B - This results because the propa- 
gator for the fictitious unphysical A "particle" enters with the wrong sign, by construction. 
However, the spirit of the regularization is to produce a well-defined effective theory, valid 
up to the scale of the fictitious PV "particle" A ^> Mb- The higher the value of A, the more 
stable is the effective theory. We are therefore led to ask, at what value of A/M B does the 
effective theory for the "neutrinos only" ansatz fail to support the idealized "neutrinos only" 
model? We answer this question with the graphs displayed in Figs. 1 and 2. If the theory 
fails at A not much larger than 2M X , then the field theory is very incomplete as written; 
more fields are needed already near the dark matter scale M x . With more fields come more 
counterterms. Consequently, more fine tuning of these counterterms is needed to maintain 
the "neutrinos only" model. It seems more likely that nature eschews this additional tedium, 
and the electromagnetic branching ratio increases at scales > A. 

In Figs. 1 and 2 we show the ratio 1Z of annihilation rates into l + l~ (through the box 
diagram) versus into vv, as a function of the cutoff scale (in units of M B ), with M B taken to 
be 10 and 20 times the W mass, respectively, and for M x taken to be 2, 5, and 10 times M w . 
In calculating the box graph, we have neglected complicating non-leading contributions, as 
they are unnecessary for the thesis of this paper. From the figures, one sees that already at 
A ~ 3 times Mb, the relative rate into the l + l~ mode is of order 1%. In fact, we may invoke 
the simple asymptotic formula in Eq. (19) to infer that the rate ratio 1Z grows as ln 2 (A/M B ), 
and that a somewhat minimal mass ordering A ~ 10 Mb ~ 10 2 Mw will necessarily yield 
1Z ~ few percent. 

Our results for the ratio 1Z of rates versus the x mass are shown in Figs. 3 and 4. We 
note that TZ is an increasing function of M x , of A/ Mb, and of Mb/Mw, as shown explicitly 
in Eq. (17). We note that the monotonic increase of the rate ratio with A is an inevitable 
consequence of PV regularization, and the increase with M x and M B derives from the 
propagator suppression of the tree-level rate, ~ (t — M B Y 2 ~ (— M 2 — M B )~ 2 - Fortunately, 
it is precisely the most physically-reasonable values for these parameters that yield the 
largest electromagnetic branching ratio. For a reasonable renormalization scheme, we would 
want A M B - Further, M B /M W must be large enough to avoid heretofore undetected new 
physics [10]. 

It is also interesting to determine the constraints which apply to our parameters if x is 
a thermal relic. In this case, Mb and M x are constrained by the requirement that the x 
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FIG. 1: The ratio of the rate for XX from the box diagram to the rate for XX —* vv at tree 

level, as a function of A /Mb, where A is the cutoff scale of the effective theory and Mb is the mass 
of the boson that mediates the x annihilation into neutrinos. Here, Mb/M\y = 10, and M x /My/ 
has the values 2, 5, and 10. 
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FIG. 3: The ratio of the rate for XX —* from the box diagram to the rate for \X vv at 

tree level, as a function of the % mass, for Mb/Mw = 10 and A/Mb = 2, 5, and 10. Mb is the 
mass of the boson that mediates the x annihilation into neutrinos, and A is the cut-off scale of the 
effective theory. 
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FIG. 4: As Fig. 3, but with M B /M W = 20. 
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abundance account for the currently observed dark matter: the value fl x h 2 = 0.1 requires 
an annihilation rate of approximately (va) ~ 3 x 10" 26 cm 3 sec -1 for thermally-produced 
dark matter [1, 11]. In the model presented here, in the non-relativistic limit s ~ (2M X ) 2 
and large M\ ^> s limit, we have 



(vv)xBr(xx^™)= (!?) (t^) (2) 



Ml 



4tt J \M% 

for the pure L = partial wave s-channel exchange of the B, and ~ \ times this for t- 
channel B exchange (where an L = 1 contribution to the rate is suppressed by the square 
of the x velocity); Br(xx ~~ is the annihilation branching ratio to vv. Substituting the 
desired rate for (va) into this equation, and taking Br(xx ~^ — 1) we obtain the result 

2 

M x = 70 TeV, (3) 



4tt\ fMi 

7b) \m : 



Mb (4) 



or 

M x " V \4ir J V M x 
If we consider only M x > 100 GeV, and restrict the coupling to the perturbative regime 
9b — then we obtain the bound [29] 

K-^<26, (5) 

where the lower bound comes from the additional requirement that M x < Mb so that the x 
dark matter is not de-stabilized by the simple decay chain x ~^ B + v [30] . Thus, the mass 
of the boson B cannot be wildly larger than that of the dark matter particle. (When the 
box-diagram annihilation rate into l + l~ becomes significant, the value of Br(xX ~^ < 1 
must be included in the algebra.) Of course, if the dark matter is not a thermal relic, 
then these constraints no longer apply. None of the other results in this paper rely on the 
assumption that x is a thermal relic, and of course, the constraints on (va) that we examine 
later in the paper are irrelevant if (va) is fixed to give the thermal relic abundance. 

We have presented figures only for the case where M x > My/. However, the calculation 
in Appendix (I) is valid for the case M x < M w as well. In this case, we find that the 
branching into l + l~ is infinitesimal. However, as we have not included the induced Zxx 
vertex shown in Fig. (c) in the calculation of the annihilation rate into v&, we cannot offer 
a reliable estimate of the rate if M x < M w . A calculation of the contribution of the vertex 
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diagram to the total annihilation rate is irrelevant to conclusions of the present paper, but 
would be needed to extend the framework presented here to a more complete "neutrinos 
only" annihilation model for M x < M w . 

Returning to the case M x > Mw, our results indicate that a model for dark matter 
annihilation that produces only neutrinos at tree level, when sensibly embedded into a field 
theory, will produce a significant (> 1%) branching into . It is reasonable to ask whether 
such a branching ratio is large enough to provide limits competitive with the neutrino limits 
discussed in Refs. [3, 4, 6]. The neutrino annihilation limits in Ref. [3] are for cosmic 
dark matter annihilations, while the (tighter) limits from Ref. [4] are for annihilations in 
the halo of the galaxy; it is the latter which we will consider here. (Palomares-Ruiz and 
Pascoli [5] consider neutrino constraints for dark matter masses below 100 MeV, outside of 
the mass range for which our results are interesting). Ref. [6] also provides general limits 
on annihilation into high-energy photons. These limits are not directly applicable to our 
calculation, since we consider branching into charged leptons rather than photons. High- 
energy leptons will convert their energy into photons primarily through inverse Compton 
scattering off of background photons, or, if a magnetic field is present, through synchrotron 
radiation. For our purposes, it is sufficient to consider only the latter effect, within the 
Galaxy, where the magnetic field is reasonably well known. For discussions of the former 
effect, see, e.g., Refs. [12-16]. (See also Ref. [17] for one of the earliest discussions of these 
phenomena) . 

Synchroton radiation from the products of dark matter annihilation has been investigated 
in some detail [12-14, 18-23]. More recently, it has been suggested that the microwave 
"haze" near the center of the Galaxy, observed by WMAP, could be synchrotron radiation 
from such annihilations [24-26]. Finally, Hooper [27] has used the WMAP observations to 
bound annihilation-produced synchrotron radiation. This latter result can be combined with 
our predicted branching ratio into charged leptons and compared with the corresponding 
neutrino bounds. 

Our results allow branching into e + e~, or t + t~, depending on the neutrinos that 

couple in £ x bu- Synchrotron radiation is produced directly by e + e~, and indirectly through 
the decay-produced electrons from ji + /i~ and t + t~. In the latter two cases, some of the 
decay energy goes into neutrinos and does not contribute to the synchrotron radiation. Thus, 
the e + e~ annihilation channel yields the strongest synchrotron signal, followed by fi + fi~ and 
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then t + t~ [25]. Since we make no prediction for the specific type of neutrinos produced 
by the dark matter annihilations, we can derive the most conservative limits by considering 
annihilation into t + t~ only. 

Hooper [27] gives limits on the cross-section for dark matter annihilation into t + t~ using 
the requirement that the resulting synchrotron radiation not exceed the WMAP microwave 
observations near the center of the Galaxy. He presents two sets of limits, the first cor- 
responding to an assumed Navarro- Frenk- White (NFW) profile [28], and the second to an 
assumed uniform distribution of the dark matter inside of the solar circle. Similarly, Yuksel 
et al. [4] give three different sets of limits based on the angular scale over which the neu- 
trino emission is measured, and they provide appropriate modification factors for each limit 
assuming any of three density profiles. While there is no exact correspondence between any 
of the limits provided in Ref. [27] and those given in Ref. [4], the closest correspondence is 
between the NFW sychrotron limit in Ref. [27] and the "Halo Angular" limit in Ref. [4], 
for the parameters appropriate to the NFW profile. Comparing these two sets of limits for 
annihilation purely into t + t~, with M x = 100 GeV, we find that the upper bound on (va) 
from vv production is roughly 1000 times the corresponding limit from r-produced syn- 
chrotron radiation. This ratio decreases to about 30 at M x = 1 TeV. Thus, the synchrotron 
bounds from the production of charged leptons are tighter than the neutrino bounds at 
M x ~ 100 GeV if the branching ratio to charged leptons is > 0.1%. For M x ~ 1 TeV, this 
branching ratio must be > 3% in order for the synchrotron bounds to be tighter than the 
neutrino bounds. This is the mass range (100 GeV < M x < 1 TeV) of greatest interest from 
the point of view of WIMP dark matter. (Strictly speaking, the neutrinos from r decay 
should be included in the total neutrino signal when comparing to the neutrino limits in 
Ref. [4], but this is a small effect as long as the branching ratio into charged leptons is 1). 

Our results do not undercut the neutrino limits proposed in Refs. [3-5]; these certainly 
remain valid limits on the annihilation rate for dark matter. Further, it is clear that both 
the synchrotron limits derived in Ref. [27] and the neutrino limits in Ref. [4] depend on the 
assumed model for the distribution of dark matter in the Galaxy. However, our purpose 
in this paper is not to determine the exact ratio between these two sets of limits (which 
may change with improved observations or more detailed calculations in any case) but 
simply to point out that the neutrino and charged lepton/photon sensitivities are likely 
to be competitive with each other over the WIMP mass range of greatest interest. For dark 
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matter masses below the W mass, we can offer no useful new limits. Our calculations in this 
case allow annihilation into neutrinos only, with only an infinitesimal branching to other 
particles. 

Although our results might seem to be limited to the specific model we have examined, 
they are rather general. Whenever neutrino-antineutrino pairs are produced in the final 
state, they can always be converted into l + l~ pairs through Z or W exchange as in diagrams 
(a) to (d) above. This is an inevitable consequence of the Standard Model; no new physics 
is required. The only assumption concerns the creation of a specific interaction that couples 
dark matter to vv pairs at tree level. Once such a model is in place, however, the induced 
coupling to l + l~ is inevitable and readily calculable. 
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I. APPENDIX: CALCULATION OF BOX DIAGRAM (IN UNITARY GAUGE) 
AND ELECTROMAGNETIC RATE 

Here we present the amplitude for the box diagram of Figure (d), and a calculation of 
the divergent part of this amplitude. Then we construct the ratio of rates, 



The subscripted qualifier "DivPart" is a reminder that only the leading (divergent) contri- 
bution to the electromagnetic rate will be included here. 
In unitary gauge, the IF-propagator is 



A V(j (xx -»• v v ))_ 



DivPart 



(6) 



k 2 -M^ + iM w T w ' 
and the matrix element for the amplitude of the box (d) is given by 



(7) 
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x [u(q)T^ + £ + M v )Y B u(p)\ [v(p)r B (i - t + M p )T» w v(q)] , 

where the momenta assignments to the particles are x(p)i x(p)> l~(o), ^ + (<?)> W(k), which 
in turn determine the further assignments B(k — p + q), v[k + q), and P(q — k). Here 
= (-^)7 A ' ( - 1 ~ 75 ' > is the usual charged-current electroweak vertex with the SU{2) coupling 
constant g — e/ sind m , p is the four-momentum of the dark matter, q is the four-momentum 
of the electron, and the internal loop four-momentum of the virtual lU-boson is given by k. 
In addition, in the numerator is the coupling times Lorentz structure assigned to a B\v 
vertex, not to be identified with the B width in the denominator. The value of the coupling 
we do not require, as the coupling cancels out when divided by the amplitude of the tree 
diagram; accordingly, we let simplicity be our guide and choose a scalar Lorentz structure 
for the B field. We will neglect the neutrino and charged lepton masses, and define M x to 
be the dark matter mass. 

We introduce four Feynman parameters 6 (one per internal line) and find for the denom- 
inator D: 

^ = (4 - 1)! T d& f di 2 f df s f dU 6(1 - [6 + 6 + 6 + &]) (9) 
u Jo Jo Jo Jo 

x [&((* - P + qf - Ml) + Uk 2 - Mb) + Uk + qf + W ~ if] 

Here we have dropped the B and W widths from the denominator, as they play no important 
role in the present calculation. Upon rotating k to Euclidean space, we find that the divergent 
part of the amplitude, coming from the longitudinal mode of the W propagator, is 

M ^ = ^rj^iJ d Z {k 2 E + A 2 y x EfoX 1 + 7 5 )rWp)] [v(P)^b(i - TbMg)] ■ (io) 

Here we have defined the functional 

/ di = 3! f di, f d& f de 3 f dU S (1 - [6 + 6 + 6 + 6]) (11) 
J Jo Jo Jo Jo 

/ dCi / d6 / 

LJo Jo Jo J &=i-£i-&-& 



= 6 



and the parametrized mass-squared 

A| = 6 Mi + 6 - (6 - ^i 2 - 266 + 466) . (12) 
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In (12) we have set the invariant t = (q — p) 2 equal to the value appropriate for a non- 
relativistic XX annihilation, namely, t ~ — M 2 . 

Performing the integral over k B will result in a logarithmic divergence, and therefore we 
introduce a Pauli-Villars regulator for the B propagator 

i i i 



(13) 



(k-p + q) 2 -M 2 (k-p + q) 2 -M 2 (k - p + q) 2 - A 2 
As is well-known, the Pauli-Villars regularization preserves local and global symmetries of 
the interaction. We then perform the integral over k B , and find that the divergent amplitude 
after regularization becomes 



M box = 



V^M 2 W 



[«(g)(l + 75)r B «(p)] [u(p)r B (l-75 



Hq)} I 



A 2 

^A 



A 2 



(14) 



where 



A\ = A| (Mi -> A 2 ) = 6 A 2 + 6 - (6 - £ - 2^6 + M 2 



(15) 



Due to the assumed scalar nature of the B particle, the 75's in Eq. (14) may be omitted. 
The tree-level amplitude for XX ~^ v ^ is 



iM 



tree 



[u{q)V B u(p)\ 



t-M 2 



[v(p)T B v(q)} . 



(16) 



Dividing the box amplitude by the tree-level amplitude (with t ~ — M 2 , again) and squaring, 
we arrive at the result 



K 



M 2 +M 2 
128 • 8% ■ Myy 



A 2 



A 2 

B 



For the running SU (2) coupling we use the value applicable at the weak scale: 

- 2 16tt 
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sin 2 0„ 



4e = 167TQ; f« 



128 



(17) 



(18) 



Notice that the spin-parity assignment of the _B-meson appears only via the Dirac structure 
T B , which cancels out of the ratio 1Z. Thus, our results to follow are independent of the 
5-meson's spin and parity. 

There are some interesting, exactly calculable limiting cases. For M B ^> M 2 , M^, one 
has A\/A B ~ A 2 /M B: independent of J d£, and so Eq. (17) becomes to lowest non-vanishing 
order, simply 



K — ► 0.97 x 1(T 7 



M 



B 



M 



w 



In 



A^ 
M| 



for Ml > M 2 Mw . 



(19) 
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We note that Eq. (5) allows this mass ordering, but does not require it. 

Another interesting calculable limit is M 2 <C M|,A 2 ,M^. Neglecting M 2 in A 2 B and 
A^, one finds that to lowest non- vanishing order, 



n 



Ml 



128 • 8tt • M 2 , 



In 



A^_ 

Ml 



+ 



0-4) 



for M 2 <Ml,A 2 ,Mj 2 F . 

(20) 



If in addition to M 2 < Mj,,A 2 ,M^, one includes M|,A 2 < M^, i.e., a low-mass model 



with a low-mass cutoff, then there results 



K -> 0.97x10 



-7 



M 



Mi 



A 5 



In 



A 2 



Ml 
M 2 W 



In 



M 2 



1 2 



, for M 2 < M|, A 2 < M, 
(21) 
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Another limit, for a low- mass model with a high-mass cutoff, is 

2 



TZ^0.97x 1Q- 7 



M 2 + M 2N :? 



In 



MA 



valid for M 2 M| < A 2 , 



(22) 

The bracketed quantity in (22) is a monotonically increasing function of A/M w , equal to 
zero at A = (as it must at A = M B ), to 1 at A = M W) and growing as ln(A 2 /M^) at 
A 2 3> Myy- It is clear that with the ordering Ml <C M^, i.e., when the tree graph has a 
light-mass (Mb) propagator while the box graph has a heavy-mass (My/) propagator, then 
the resulting electromagnetic branching fraction is negligible for any value of the effective 
field thoery cutoff A, even up to the Planck mass. 



II. APPENDIX: CALCULATION OF BOX DIAGRAM IN GAUGE 

Given that the unitary gauge admits only the three physical W spin states into the W 
propagator, this gauge choice is guaranteed to produce diagrams whose cuttings respect 
the Cutkowsky rules; other gauges will not. Nevertheless, the "renormalizable gauges", 
rather than the unitary gauge, are often invoked because their resulting propagators offer 
a more benign high-energy behavior. Of course, when a gauge invariant set of graphs is 
summed to calculate an S- matrix element, all gauges must give the same gauge-invariant 
(£- independent) answer. However, the box graph of interest here is not itself gauge invariant, 
the reason being that the dark sector breaks electroweak gauge invariance by treating the 
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neutrino differently from its SU(2) charged lepton partner. Consequently, we expect the 
box graph to yield a gauge-dependent answer. We have dealt with this situation in the main 
text by arguing that the unitary gauge is singled out since it, and it alone, includes only 
physical W states in its Feynman rules. Nevertheless, it is interesting and instructive to 
calculate the box graph in a general gauge. We do that here. 

The unitary gauge W^-propagator, given in Eq. (7), is the £ 
i?£-gauge VF-propagator (we neglect the finite VF-width) 

0' 



oo limit of the general 



and the accompaning unphysical Goldstone boson propagator is 

+i 



kfjky (1 



(23) 



D GB (0 



(24) 



In contrast to the calculation in the unitary gauge where the leading high-energy behavior 
was given by a logarithmic divergence, here, for any finite 6 the leading high energy con- 
tribution to the box graph is finite and therefore calculable. However, the result will be 
gauge (£) dependent. To facilitate the calculation, we note that (i) the Goldstone boson 
contribution to the amplitude does not contribute to the leading high-energy behavior, and 
so will be neglected, and (ii) the R$ gauge W^-propagator (Eq. 23) can be rewritten as a sum 
of the ^-independent unitary gauge propagator and a ^-dependent correction: 



K(0 



+ 



k 2 - £M 2 



(25) 



w 



As far as the high-energy behavior of the amplitude is concerned, the contribution from the 
unitary gauge W^-propagator again yields A^ in Eq. (12), which we rename A 2 ^ to emphasize 
the role of the physical W^-propagator in this contribution. The additional contribution from 
the "wrong sign" propagator with a pole at the Goldstone boson mass-squared yields 



A 2 GB = £i M B + 6 (£M^) - (6 - H - 266 + 466) M 2 



(26) 



We then have 



K 



M 2 +M 2 
128 • 8tt • M 2 , 



The bracketed expression here, f d£ In -^f 2 - , is bounded by In 
the fictitious Goldstone boson mass-squared has come to play the role in the ii^-gauge 



d£ln 



A 2 

GB 



A 2 
w 



^7 



(27) 

In |£|. We see that 
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calculation that the fictitious PV mass-squared A 2 played in the unitary gauge calculation. 
As a consistency check, we see that both £ — > oo and A 2 — > oo return the infinite unitary 
gauge contribution without PV regularization. 

This result explicitly demonstrates the merit of using the unitary gauge and a Pauli- 
Villars regularization scheme for the divergent graph. The particle degrees of freedom in 
the unitary gauge are all physical, and the PV cutoff A has physical meaning. A is the scale 
up to which the model is applicable, and beyond which the ultraviolet completion of the 
theory cannot be ignored. In contrast, the cutoff provided by the Goldstone mass in the 
renormalizable gauges has no physical interpretation. 
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